Roundabout relaxation: collective excitation requires a detour to equilibrium 
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Relaxation to equilibrium after strong and collective excitation is studied, by using a Hamiltonian 
dynamical system of one dimensional XY model. After an excitation of a domain of K elements, the 
excitation is concentrated to fewer elements, which are made farther away from equilibrium, and 
the excitation intensity increases logarithmically with K. Equilibrium is reached only after taking 
this "roundabout" route, with the time for relaxation diverging asymptotically as K 1 with 7 ~ 4.2. 

PACS numbers: 05.70.Ln, 05.45.-a, 87.10. +c 



43 

o 



CO 



I 

o 
o 



> 
o 

o 
o 



-a 
c 

o 
o 



X 



Relaxation to equilibrium has been one of the most 
important topics in non-equilibrium phenomena. The 
relaxation dynamics after weak excitation has been thor- 
oughly investigated; it is represented as a superposition 
of the dynamics of excited modes, each of which depends 
on external parameters such as temperature. When sev- 
eral modes are strongly and collectively excited, on the 
other hand, interaction among the excited modes is not 
negligible, which may cause non-trivial dynamic behav- 
iors that are not simply represented by the superposition 
of the excited modes. In particular, a set of the excited 
modes may form a partial system, and give an internal 
state that cannot be determined solely with an external 
conditions. The relaxation will depend on the internal 
state, and, in turn, the internal state will dynamically 
change with the relaxation; such an interplay between the 
internal state and the relaxation is expected to be seen 
in systems far from equilibrium in general. We intend 
to search for some non-trivial and universal relaxation 
phenomena therein. 

Previously the authors have reported such a novel 
relaxation phenomenon, using a Hamiltonian dynami- 
cal system of XY model with mean field coupling [1]; 
when a part of the system is highly excited, the relax- 
ation progresses intermittently through bottlenecks by 
self-organizing a critical state for the partial excited sys- 
tem. Now it is interesting to study the relaxation in 
a corresponding lattice system. In particular we study 
a one dimensional (ID) XY model, which has no phase 
transition and accordingly no critical state. In spite of its 
absence, we find a rather remarkable relaxation process, 
i.e. "roundabout" relaxation that a partially excited sys- 
tem reaches equilibrium only after it once goes farther 
away from equilibrium. In the present Letter, we report 
this discovery and analyze its mechanism, by emphasiz- 
ing the divergence of the relaxation time to equilibrium 
with the number of excited elements. 

The Hamiltonian that we study is [2] 
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where 6i G [— it, 7t), with periodic boundary condition 
0i = 6i+N- We consider the case of J > (ferromag- 
netic), setting J = 1 without loss of generality. The N 
pcndula, placed on the ID lattice, are coupled by the 
nearest neighbor interaction. Each pendulum has two 
types of motion; rotation at a higher and libration at a 
lower energy. 

Note that when a single element is highly excited, its 
relaxation is slow, with the relaxation time increasing ex- 
ponentially with the excitation momentum, as in the case 
with mean field coupling [3]. This is simply because the 
effective interaction between the excited element and the 
others rapidly decreases as the excitation momentum in- 
creases. The exponential dependence as such is discussed 
in terms of Boltzmann- Jeans conjecture (BJC) [4]. 

Here we study the relaxation when a domain of the sys- 
tem consisting of quite a few elements is highly excited, 
while the whole system is much larger [1]. To be specific, 
the system is prepared to be in equilibrium, with an en- 
ergy density U and total momentum zero. At t = 0, 
K elements are simultaneously excited with the same 
momentum Pk- The momentum profile then becomes 
almost rectangular, with two domains separated by two 
interfaces. The number of the excited elements are much 
smaller than that of the total system, K/N C 1; we keep 
K/N = 0.1 in the following numerical simulation. 

We first numerically observe the relaxation process. 
The element at the interfaces of the excited domain is 
intermittently absorbed into the not-excited domain, los- 
ing its excited energy. The population of the excited el- 
ements accordingly decreases one by one from the 
initial K toward zero (equilibrium). If the excitation is 
strong enough, the excited domain receives some posi- 
tive momentum from the escaped element, to increase 
its center-of-mass momentum (CMM); the process is re- 
peated as the relaxation progresses (FIG. 1). The excited 
elements then rotate faster and faster, and thereby go far- 
ther away from equilibrium. The increase of the CMM 
implies more time for the next escape, as in the argu- 
ment for BJC, and requires much longer time to reach 
equilibrium. 

Thus the highly excited state is self-sustained, and the 
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FIG. 1: Two snapshots of the momentum profile p i , with 7 
as the time average over a period of 1000, at t=0 and 25000 
(corresponding to the arrows to the inset). In the inset, the 
corresponding time series of the momenta of the excited ele- 
ments are overlaid, where abscissa axis is logarithmic scale in 
time. N = 100, K = 10, P K = 6, U = 0.3. 
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FIG. 2: The change of the center-of-mass momentum 
(x, solid line) and the effective temperature AT' E ' (+, dotted 
line), which are time averaged over a period of 100, against the 
(initial) effective temperature Tr. N = 300, K = 30, Pk = 8. 
Obtained from the numerical calculations of 25 samples for 
each T K . 



relaxation to equilibrium must take this "roundabout" 
route. This process is not observed near equilibrium, and 
requires sufficient excitation of a domain of elements. 

Both of the two domains are almost in equilibrium in 
the inertial frame of their own CMMs, since their ele- 
ments interact mostly among themselves. The interac- 
tion at the interface, on the other hand, is much weaker, 
because the difference between each CMM is much larger 
and the interaction decreases rapidly with it. Since both 
domains are composed of a sufficient number of elements, 
each of which is regarded as a thermodynamic system, 
and accordingly the whole system as the two thermody- 
namic systems weakly coupled through the interface. 

Being large enough, the non-excited system plays the 
role of a heat bath, whose thermodynamic state is kept 
almost constant. The state of the excited domain, on the 
other hand, changes with relaxation, in the time scale 
much slower than that for equilibration within the do- 
main but much faster than the whole relaxation to the 
equilibrium. Thus it is possible to define an internal 
state for the excited system as its thermodynamic state. 
To quantify the internal state, we introduce the effective 
temperature of the excited domain, 

TiE) = ^£%- p<B5 ) a . ( 2 ) 

i 

in the inertial frame with the following CMM of the ex- 
cited domain, 

i 

where denotes summation over only the excited el- 
ements. 



To study the dynamics of the internal state, we calcu- 
late the change of the CMM, AP (B) , and of the effective 
temperature, of the excited domain, due to the 

first absorption of an interfacial element into the non- 
excited (bath) part (K — » K — 1). Their dependence on 

Tjc d = T (E) |t=o is plotted in FIG. 2, which is computed 
by preparing the excited state with from a thermal 
equilibrium state with the CMM Pk- The figure shows 
that decreases toward from the positive side, 

indicating that increases to approach a certain tem- 
perature. On the other hand, AP® approaches a certain 
positive value. Except for the K dependence of the limit 
value of AP^ , the other properties above are indepen- 
dent of any other parameter including U and sufficiently 
large Pk- By noting that single elements are absorbed 
successively in the relaxation course, the above result in- 
dicates that the effective temperature increases and ap- 
proaches a high value, and there the CMM increases by 
some constant value per each absorption. The CMM of 
the excited domain thus continues to increase, until the 
last element is absorbed, while a high temperature at the 
excited domain is sustained. 

Now we study the increase of the relaxation time with 
the number of the excited elements. Numerically it is not 
so easy to trace all the relaxation course, since it requires 
a huge time. Instead we give an analytical estimate for 
it, by focusing on the dynamics of the macroscopic quan- 
tities of the excited part. To make notations simple in 
the following analysis, subscripts is omitted, and let 
■fc be macroscopic quantities when the number of the re- 
maining excited elements is k. 

Firstly we study how the change of the CMM Pk-i—Pk 
depends on the population of the excited elements k. 
Pk-i — Pk with its dependence on K, obtained numer- 
ically, is plotted in FIG. 3. We have also confirmed that 
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FIG. 3: Pk-i - Pk versus K. P K = 8,T K « 2.4. Obtained 
from the numerical calculations of 25 samples for each K. 



FIG. 4: r versus P K . N = 300, if = 30, T K « 2.4. Obtained 
from the numerical calculations of 25 samples for each Pk- 



the relationship is almost not affected by Pk and Tk, as 
long as they are large enough. Thus we obtain, 



P k - 1 -P k =a/(k-l). 



(4) 



Here a is a constant with the dimension of momentum, 
which we numerically found a ~ 3. We examine a rough 
estimation of the value, by considering the dynamics of 
only three interfacial elements: an escaping element with 
two neighboring ones. Furthermore the motion of each 
neighboring element is approximated as a constant rota- 
tion with the momentum and SI. Then the effective 
Hamiltonian of the interfacial dynamics can be given as, 



Hi 



p 2 /2 + J[l - cos 6>] + J[l - cos(6» - Sit)]. (5) 



For |f2| > Slthj. w 6, two motions withing the region 
around p ~ and p ~ SI are separated by KAM tori, 
while for |f2| < Sl^r, an orbit can cross between the 
two regions, due to the collapse of the last KAM torus 
through the so-called "resonance overlap [5]." In other 
wards, the interface element can jump over from p ~ SI 
to p ~ 0, only if |f2| is less than Slthr- Coming back 
to the original problem, this is nothing but the escape 
of the interfacial element from the excited domain. The 
loss of the momentum Sl t hr thereby should be compen- 
sated by the the momentum gain of the neighboring ele- 
ments, each of which is given approximately Sl t f lr /2 from 
the symmetry. Hence the excited domain gains momen- 
tum about Sl t h r /2, which is distributed to the k — 1 el- 
ements therein, leading to the increase of the CMM by 
fW/2/(fc-l)w3/(fc-l). 

Secondly we compute the relaxation time r when the 
first one of the excited elements loses its energy and is 
absorbed into the non-excited part (K — ► K — 1), whose 
dependence on Pk is plotted in FIG. 4. This exponential 
dependence is just the form of the BJC mentioned above. 
This relationship again is almost not affected by K and 
Tk [6], as long as they are large enough. Thus we obtain, 



r(P) = CeP p . 



(6) 



The constant ft has the dimension of inverse momentum, 
which is the order of the inverse of the momentum for the 
separatrix motion of a pendulum. Numerically we have 
obtained (3 « 1.4 and C « 0.1. 

On the basis of the two short-time properties (4) and 
(6), we can estimate the long-term relaxation. First 
we consider the moment of the last excited element P\. 
When the population falls to k from initial K, the gain 
of CMM is, 
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= a(S K -i - Sk-i), 
by summing up (4), where 



5 * = Et 



(7) 



(8) 



which increases asymptotically as log if with K —> oo. 
The momentum of the last excited element P\ is then 
evaluated as, 



Pi -Pi 



K 



aS 



K-l 



log if. 



(9) 



Hence the CMM diverges asymptotically as log if. 

Finally we estimate the total relaxation time to the 
equilibrium. From (6) and (7), the total time up to the 
relaxation from k to k — 1 is given by, 



t k = Y^r{P j )=t K e^-^e-^ 

j=k j=k 



(10) 



def 



where 7 = a/3 is a dimcnsionless number. Then the 
relaxation time to equilibrium is estimated as, 



Ten =h= t K e lSK - 



K-l 
k=l 



(11) 
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Recalling log(fc + 1) < Sk < logfc + 1, for k > 1, we get 
e-Tfc" 7 < < (k + 1)~ 7 . (12) 

Combining (11) and (12) yields 

K-l AT-l 



1+ V ^<l£ie^ s — <1+V 



T, 
t K 



1 



As X -» 00, e^ 5 * -> e 7i°g-ff = #7, and 

if 7 > 1 
if 7 = 1 
if < 7 < 1, 



(13) 




(14) 



where £(•) is Riemann zeta. Hence the asymptotic form 
of the relaxation time to equilibrium r eq as K — ► 00 is 



K"< if 7 > 1 

^log^ if 7=1 
if if < 7 < 1. 



(15) 



Note that, if each excited elements relaxes indepen- 
dently, r eq — Ktx- Thus 7 = 1 is the lower bound of 
the divergence attributed to the cooperative effect of the 
excited elements. Indeed, in the present model, we have 
a « 3 and (3 « 1.4, and we get 7 « 4.2. Hence the 
relaxation time satisfies r eq — > if 42 , showing the rapid 
divergence with the number of excited elements. 

To check the validity of the above analysis, we have 
computed the increase of CMM and the relaxation time 
up to the fc-th element, Pk and tk (FIG. 5- (a) and (b), re- 
spectively). Even though it is hard to follow all the relax- 
ation course numerically, the numerical results (crosses) 
plotted against k/K rather well agree with the theoreti- 
cal estimates (7) and (10) (solid curves). 

In summary, we have discovered that a class of one- 
dimensional lattice system must take a "roundabout re- 
laxation route," once highly excited. After an excitation 
of a domain of the system, the momentum of the do- 
main increases logarithmically, to go farther away from 
equilibrium, before reaching equilibrium. The relaxation 
time accordingly diverges as a power K 1 against the in- 
crease of the excited elements K . Hence the relaxation 
with collective excitation has a rather peculiar form, and 
is rather different from that near equilibrium. Note that 
the present choice of one-dimensional XY (pendulum) 
model is not so special; actually preliminary results in 
two-dimensional case also suggest similar roundabout re- 
laxation. In physics there should be a variety of examples 
described by coupled pendula on a lattice, and our results 
will be relevant to slow relaxation in various real non- 
equilibrium systems. It is also notable that in a closed 
chemical reaction-diffusion system, a self-sustained ex- 
cited state is formed as a transient dissipative struc- 
ture [7], which leads to hindrance of relaxation to equi- 
librium. 
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FIG. 5: Temporal change of (a) Pk and (b) tk, plotted against 
k/K, which decreases through the relaxation. Numerical re- 
sults from 25 samples of direct relaxation simulations are plot- 
ted (x), while the solid curves are the analytical estimation 
from (7) and (10), respectively. N = 1000, K = 100, = 
8,T K »4.2. 
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